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Recent work on the mechanism of polymer crystallization has led to a proposal for the mechanism 
of thickness selection which differs from those proposed by the surface nucleation theory of Lauritzen 
and Hoffman and the entropic barrier model of Sadler and Gilmer. This has motivated us to 
reexamine the model used by Sadler and Gilmer. We again find a fixed-point attractor which 
describes the dynamical convergence of the crystal thickness to a value just larger than the minimum 
stable thickness, £ m in- This convergence arises from the combined effect of two constraints on the 
length of stems in a layer: it is unfavourable for a stem to be shorter than l m i n and for a stem 
to overhang the edge of the previous layer. The relationship between this new mechanism and 
the explanation given by Sadler and Gilmer in terms of an entropic barrier is discussed. We also 
examine the behaviour of the Sadler-Gilmer model when an energetic contribution from chain folds 
is included. 
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I. INTRODUCTION 

On crystallization from solution and the melt many 
polymers form lamellae where the polymer chain tra- 
verses the thin dimension of the crjastal many times fold- 
ing back on itself at each surfaceirtl Although lamellar 
crystals were first observed over forty years ago their 
physical origin is still controversial. It is agreed that the 
kinetics of crystallization are crucial since extended-chain 
crystals are thermodynamically more stable than lamel- 
lae. However, the explanations for the dependence of 
the lamellar thickness on temperature offered by the tw| 
dominant theoretical approaches appear irreconcilable 
The lamellar thickness is always slightly greater than 
^min, the minimum thickness for which the crystal is 
thermodynamically more stable than the melt; Z m ; n is 
approximately inversely proportional to the degree of 
supercooling £1 

The first theory, which was formulated by Lauritzen 
and Hoffman SQoii.after the initial discovery of the chain- 
folded crystals ,ETEi invokes surface nucleation of a new 
layer on the thin side faces of the lamellae as the key 
process. It assumes that there is an ensemble of crystals 
of different thickness, each of which grows with constant 
thickness. The crystals which grow most rapidly domi- 
nate this ensemble, and so the average value of the thick- 
ness in the ensemble, which is equated with the observed 
thickness, is close to the thickness for which the crystals 
have the maximum growth rate. The growth rates are de- 
rived by assuming that a new crystalline layer grows by 
the deposition of a succession of stems (straight portions 
of the polymer that traverse the crystal once) along the 
growth face. The two main factors that determine the 
growth rate are the thermodynamic driving force and 
the free energy barrier to deposition of the first stem in 



a layer. The former only favours crystallization when 
the thickness is greater than l m [ n . The latter increases 
with the thickness of the crystal because of the free ener- 
getic cost of creating the two new lateral surfaces on ei- 
ther side of the stem and makes crystallization of thicker 
crystals increasingly slow. Therefore, the growth rate 
passes through a maximum at an intermediate value of 
the thickness which is slightly greater than Z m i n . 

The second approach, which was developed by Sadler 
and Gilmer and has been termed the entropic barrier 
model, is based uooa, the interpretation of kinetic Mojatej 
Carlo simulationsfali3 and rate-theory calculationsllirMI 
of a simplified model of polymer crystal growth. The 
model has since been, applied to the crystallization of 
long-chain paraffins,Lj copolymer crystallizatioixJ and 
non-isothermal crystallization.li3 As with the surface nu- 
cleation approach, the observed thickness is suggested to 
result from the competition between a driving force and a 
free energy barrier contribution to the growth rate. How- 
ever, a different cause for the free energy barrier is postu- 
lated. As the polymer surface in the model can be rough, 
it is concluded that the details of surface nucleation of 
new layers are not important. Instead, the outer layer of 
the crystal is found to be thinner than in the bulk; this 
rounded crystal profile prevents further crystallization.L2l 
Therefore, growth of a new layer can only begin once a 
fluctuation occurs to an entropically unlikely configura- 
tion in which the crystal profile is 'squared-off '. As this 
fluctuation becomes more unlikely with increasing crystal 
thickness, the entropic barrier to crystallization increases 
with thickness. 

Recently, we performed a numerical-study of a sim- 
ple model of polymer crystallization.EZTtj Our results 
suggest that some of the assumptions of the Lauritzen- 
Hoffman (LH) theory do not hold and led us to propose 
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a new picture of the mechanism by which the thickness 
of polymer crystals is determined. Firstly, we examined 
the free energy profile for the crystallization of a polymer 
on a surface. This profile differed from that assumed by 
the LH theory because the initial nucleus was not a sin- 
gle stem but two stems connected by a fold which grew 
simultaneouslyO 

Secondly, we performed kinetic Monte Carlo simula- 
tions on a model where, similar to the LH approach, 
new crystalline layers grow by the successive deposition 
of stems across the growth face, but where some of the 
constraints present in the LH theory are removed — stems 
could be of any length and grew by the deposition of 
individual polymer units. We, therefore, refer to this 
model as the unconstrained Lauritzcn-Hoffman (ULH) 
model. Attempts had previously been made to study sim- 
ilar models but at a time when computational resources 
were limited, so that onijj [approximate and restricted 
treatments were possible.l23c-3 

These kinetic Monte Carlo simulations confirmed that 
the initial nucleus was not a single stem. We also found 
that the average stem length in a layer is not determined 
by the properties of the initial nucleus, and that the 
thickness of a layer iSp-jiugeneral, not the same as that 
of the previous layer I 1S H 19 I Furthermore, during growth 
lamellar crystals select the value of the thickness, £**, for 
which growth with constant thickness can occur, and not 
the thickness for which crystals grow most rapidly if con- 
strained to grow at constant thickness. This thickness 
selection mechanism can be understood by considering 
the free energetic costs of the polymer extending beyond 
the edges of the previous crystalline layer and of a stem 
being shorter than Z m i n , which provide upper and lower 
constraints on the length of stems in a new layer. Their 
combined effect is to cause the crystal thickness to con- 
verge dynamically to the value I** as new layers are added 
to the crystal. At I**, which is slightly larger than Z m i n , 
growth with constant thickness then occurs. 

Some similarities were found between the behaviour of 
our model and that of the Sadler-Gilmer (SG) model. 
For example, at low supercoolings rounding of the crys- 
tal growth front was found to inhibit growth. However, 
the mechanism of thickness selection that we found ap- 
pears, at first sight, to be incompatible with Sadler and 
Gilmer's interpretation in terms of an entropy barrier. 
Here, we reexamine the SG model in order to determine 
whether the fixed-point attractors that we found at large 
supercooling also occur in the SG model, and in order to 
understand the relationship between the two viewpoints. 



II. METHODS 

In the SG model the growth of a polymer crystal re- 
sults from the attachment and detachment of polymer 
units at the growth face. The rules that govern the sites 
at which these processes can occur are designed to mimic 



the effects of the chain connectivity. In the original three- 
dimensional version of the model, kinetic Monte Carlo 
simulations were performed to obtain many realizations 
of the polymer crystals that result. Averages were then 
taken over_these configurations to get the properties of 
the modela Under many conditions the growth face is 
rough and the correlations between stems in the direction 
parallel to the growth face are weak. Therefore, an even 
simpler two-dimensional version of the model was devel- 
oped in which lateral correlations are neglected entirely, 
and only a slice through the polymep-fpystal perpendic- 
ular to the growth face is considered One advantage 
of this model is that it can be formulated in terms of a 
set of rate equations which can easily be solved. The be- 
haviour of this new model was found to be very similar 
to the original three-dimensional model. 




FIG. 1. A schematic picture of a two-dimensional slice 
(perpendicular to the growth face) through a lamellar poly- 
mer crystal which forms the basis of the two-dimensional ver- 
sion of the Sadler-Gilmer model. The three possible changes 
in configuration allowed by the model are shown (the dashed 
lines represent the outline of the possible new configura- 
tions). 

The geometry of the model is shown in Figure |]. 
Changes in configuration can only occur at the outer- 
most stem and stems behind the growth face are 'pinned' 
because of the chain connectivity. There are three ways 
that a polymer unit can be added to or removed from the 
crystal: (I) The outermost stem can increase in length 
upwards. (2) A new stem can be initiated at the base of 
the previous stem. (3) A polymer unit can be removed 
from the top of the outermost stem. 
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The set of rate equations that describes the model 
can be most easily formulated in a frame of reference 
that is fixed at the growth face. The position of each 
stem (which is representative of a layer in the three- 
dimensional crystal) is then denoted by its distance from 
the growth face (the outermost stem is at position 1). 
The rate equations are formulated in terms of three quan- 
tities: C n (i), the probability that the stem at position n 
has length i; P n (i,j), the probability that the stem at po- 
sition n has length i and the stem at position n + 1 has 
length j] and f n (i,j), the conditional probability that the 
stem at position n + 1 has length j given that the stem 



at position n has length i. These quantities are related 
to each other by 

C n (i) = Y d Pn{i,3) (1) 

3 

and 

fn(i,j) = Pn(i,j)/C n (l). (2) 

The evolution of the system is then described as fol- 
lows: for i > 1 



^ ' = fc + Pi(* - 1,3) + k-(i + l,j)Pi(i + l,j) + k-(l,i)Px(l,i)h(i,j) - 2k+P 1 (i, j) - k'iij^iij), (3) 
and for i = 1 

^%ti =k+Ci{j)+k - {2J) p i{2J) + k - {1)1)p ^ (4) 



where k + and k~(i,j) are the rate constants for attach- 
ment or detachment of a polymer unit, respectively. 

The first three terms of Equation (^) correspond to the 
flux into state by extension of the outermost stem, 
reduction of the outermost stem and the removal of an 
outermost stem of length 1, respectively. The fourth term 
is the flux out of due to the extension of the outer- 
most stem or the initiation of a new stem, and the fifth 
term is the flux out due to a reduction of the outermost 
stem. The only difference for i — 1 is that the first term 
now corresponds to flux into state (l,j) by initiation of 
a new stem. _ 

Sadler and Gilmer found /„ to be independent of n.til 
Therefore, fa can be replaced by /i in Equations (||) and 
(Q) and so the steady-state solution (dP(i, j)/dt = 0) 
of the above c max x c max equations can easily be found 
(cmax is the size of the box in the c=direction) . We used a 
fourth-order Runge-Kutta schemcEj to integrate the rate 
equations to convergence starting from an initial config- 
uration that, for convenience, we chose to be a crystal of 
constant thickness. From Pi it is simple to calculate all 
P n by iterative application of the equations 

Cn+l(j) = X)P»(;,j) (5) 

i 

and 

Pn+i{i,3) = f{i,j)C n +i(i)- (6) 

Far enough into the crystal the properties converge, i.e. 
C n +i = C n . We denote the quantities in this limit as 
Cbuik and Pbuik- 

The rate constants, fc + and k~(i,j) are related to the 
thermodynamics of the model through 

k + /k~ = exp(-AP/fcT), (7) 



where AP is the change in free energy on addition of a 
particular polymer unit. The above equation only defines 
the relative rates and not how the the free energy change 
is apportioned between the forward and backward rate 
constants. We follow Sadler and Gilmer and choose k + 
to be constant. We use l/k + as our unit of time. 

In the model the energy of interaction between two 
adjacent crystal units is -e and the change in entropy 
on melting of the crystal is given by AS* = AH/T m = 
2e/T m , where T m is the melting temperature (of an in- 
finitely thick crystal) and AP is the change in enthalpy. 
It is assumed that AS is independent of temperature. 
We can also include the contribution of chain folds to 
the thermodynamics by associating jtiie energy of a fold, 
e f , with the initiation of a new stem.Ej Sadler and Gilmer 
ignored this contribution and so for the sake of compar- 
ison most of our results are also for e/ = 0. 

From the above considerations it follows that the rate 
constants are given by 

k-(i,j) = k + exp{2e/kT m -e/kT + e f /kT) i = l (8) 
k~(i, j) = fc+exp(2e/fcT m - 2e/kT) i<j,i^l (9) 
k~(i,j) = fc + exp(2e/fcT m - e/kT) i > j. (10) 

The first term in the exponents is due to the gain in en- 
tropy as a result of the removal of a unit from the crystal; 
the second term is due to the loss of contacts between the 
removed unit and the rest of the crystal and the third 
term (only for i = 1) is due to the reduction of the area 
of the fold surface of the crystal. 

Quantities such as the growth rate, G, and the aver- 
age thickness of the n th layer, l n , can be easily calculated 
from the steady-state solution of Pi(i,j): 

G = k + -fc-Ci(l), (11) 

where k~ is given by Equation (j^), and 
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In = }jC n (j). 



(12) 



(a) 35 



We can also estimate the minimum thermodynamically 
stable thickness, l m i n . If we assume the crystal to be of 
constant thickness: 



2eAT 



(13) 



where AT = T — T m is the supercooling. However, this 
value is a lower bound since the roughness in the fold 
surface reduces the surface free energy; the resulting sur- 
face entropy more than compensates for the increase in 
surface energy. 

As well as solving the rate equations of this model, we 
sometimes use kinetic Monte Carlo. This method is more 
appropriate when examining the evolution of the system 
towards the steady state£3 At each step in the kinetic 
Monte Carlo simulation we randomly choose a state, b, 
from the three states connected to the current state, a, 
with a probability given by 



Pab = 



k a b 



Zjfe' kab 



and update the time by an increment 

logp 



At 



J2b kab 



(14) 



(15) 



where p is a random number in the range [0,1]. 



III. RESULTS 

A. Mechanism of thickness selection 

The current model has three variables: kT m /e, T/T m 
and e//e. In this section we only discuss results for e/ = 
as we wish to reexamine the same model as that used by 
Sadler and Gilmer. Throughout the paper we chose to 
use kT m /e = 0.5. Sadler and Gilmer found that this pa- 
ramete*-,did not affect the qualitative behaviour of the 
model J13 

Figures |^ and |3| show the typical behaviour of the 
model. The crystal thickness, lbu\k, is a monotonically 
increasing function of the temperature which goes to oo 
at T m . The approximately linear character of Figure ||b 
is in agreement with the expected temperature depen- 
dence of the growth, rate: G oc exp Kg /TAT, where 
Kg is a constant .Bc3 The rounding of the edge of the 
lamellar crystal is clear from the typical crystal profile 
shown in Figure 3 a. This rounding increases with tem- 
perature (Figure iz) suggesting that it is associated with 
the greater entropy of configurations with a rounded pro- 
file. The sounding also increases with an increasing value 
of kT m /eB 
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FIG. 2. The temperature dependence of the (a) crystal 
thickness and (b) the growth rate. The lines are labelled by 
the value of e//e. 



The probability distributions of the stem length, C„, 
of course, reflect this rounding (Figure [|b). The main 
difference between C n for different values of n occurs at 
small stem lengths. All the C„ seem to have the same 
approximately exponential decay at large stem lengths. 
However, the bulk can tolerate few stems shorter than 
^min for reasons of thermodynamic stability. Near the 
growth face this thermodynamic constraint is more re- 
laxed. Moreover, at the growth face there is the kinetic 
effect that each new stem must initially be only one unit 
long, thus contributing to Ci(l). It is also interesting to 
note that Cbuik is approximately symmetrical about its 
maximum. 

The rounding has a major effect on the mechanism of 
growth. Since the majority of the short stems at the 
surface cannot be incorporated into the bulk, rounding 
inhibits growth. Therefore, a fluctuation in the outer 
layer to a thickness similar to ?buik is required before the 
growth front can advance. By this mechanism stems with 
length less than l min are selected out. From Figure ||a 
it can be seen that this selection process becomes more 
complete as one goes further into the crystal. 
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FIG. 3. (a) The thickness profile of the crystal and (b) 
the probability distribution of the stem length in layer n at 
T = 0.95T m . The lines are labelled by the value of n. (c) 
The temperature dependence of the rounding of the crystal 
profile as measured by Zi/Zbulk- The lines are labelled by the 
value of e//e, and those without data points are derived from 
Equation (|25[). 

One of the main aims of this paper is to discover 
whether the mechanism of th ickn ess determination that 
we found for the ULH modelt3'll3 also holds for the SG 
model. Central to this mechanism was a fixed-point at- 
tractor that related the thickness of a layer to the thick- 
ness of the previous layer. As layers in that model once 



grown could not thereafter change their thickness, the 
relationship between the thickness of a layer and the pre- 
vious layer was independent of the position of the layer 
in the crystal. The same is not true of the SG model 
since there is a zone near to the growth face where the 
layers have not yet achieved their final thickness (Figure 
. Therefore, it is appropriate to look for the attractor 
in the relationship between the thickness of layers in the 
bulk of the crystal. 

To do so we define f' n {i,j) as the conditional probabil- 
ity that the stem at position n has length i given that 
the stem at position n+l has length j. f' n {i,j) is simply 
related to quantities that we know: 

fUhj) = Pn(i,j)/C n+ i(j)- (16) 

From this we can calculate l n (j), which is the average 
thickness of layer n given that layer n + l has thickness 
j, using 

ln(j) = Y^fUh3)i- (17) 

i 

The actual thickness of the n th layer is simply related to 
this function by 

r„ = ^c n+1 (i)/ n (i). (is) 

j 

Figure ||a shows the function ibulk(i) at T/T m = 0.95. 
As anticipated, the figure has the form of a fixed-point 
attractor. The arrowed dotted lines show the thickness 
of the layer converging to the fixed point from above and 
below. For example, a layer in the bulk of the crystal 
which is thirty units thick will be followed by one that is 
17.8, which will be followed by one that is 14.3, . . . , until 
the thickness at the fixed point, I**, is reached. I** is 
approximately equal to the crystal thickness, /bulk- Fig- 
ure shows the crystal profile that results from starting 
kinetic Monte Carlo simulations from a crystal that is ini- 
tially 30 units thicks. The convergence of the thickness 
to I** is apparent and the resulting step on the surf a re. 
is similar to that produced by a temperature jumpE3'EHl 
Figure ^b shows that similar maps with a fixed point oc- 
cur at all temperatures. As the temperature increases 
the slope of the curves become closer to 1 and so the 
convergence to I** becomes slower. 

To understand why the function l n (J) has the form of 
a fixed-point attractor we examine probability distribu- 
tions of the stem length when the previous layer had a 
specific thickness. In Figure || we show /b u ik(*>i) l° r dif- 
ferent j. These probability distributions have a number 
of features. Firstly, the probability that the stem length 
is smaller than Z m ; n drops off exponentially. This is a 
simple consequence of the thermodynamics. Secondly, 
the probability that the stem length is larger than the 
thickness of the previous layer also drops off exponen- 
tially. This is again related to the thermodynamics. The 
absence of an interaction of the polymer with the sur- 
face makes it unfavourable for a stem to extend beyond 
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the edge of the previous layer. This feature also causes a 
build-up in probability at stem lengths at or just shorter 
than the thickness of the previous layer. Therefore, as we 
saw in our previous modelJIS'tJ there is a range of stem 



lengths between Z m i n and the thickness of the previous 
layer which are thermodynamically viable, and it is the 
combined effect of the two constraints that causes the 
thickness to converge to the fixed point. 
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FIG. 4. (a)-(c) The dependence of the thickness of a layer in the bulk of the crystal on the thickness of the previous 
layer, j, and the temperature. In (a) the fact that the mapping is a fixed-point attractor is illustrated by the dotted arrowed 
lines which show the thickness converging to the fixed point from above and below. T = 0.95T m . In (b) the lines are labelled 
by T/T m and in (c) the lines are labelled by the value of j. (d) the average crystal profile that results from an initial crystal 
that is 30 units thick at T = 0.95T m . The position of the edge of the initial crystal is at zero. 



Thirdly, when the previous layer is thick the probabil- 
ity also falls off exponentially with increasing stem length 
in the range / m i n < i < j. This decay of the probabil- 
ity is much less rapid than for stems that overhang the 
previous layer and is a kinetic effect. At each step in the 
growth of a stem there is a certain probability, p ne w, that 
a new stem will be successfully initiated, thus bringing 
the growth of the original stem to an end. It therefore 
follows that the probability will decay with an exponent 
log(l — Pncw)- It is for this reason that the thickness 
of a new layer will remain finite even if the previous is 
infinitely thick. Comparison of /bulk at different tem- 



peratures shows that this rate of decay increases with 
decreasing temperature. This is because the thermody- 
namic driving force is greater at lower temperature mak- 
ing successful initiation of a new stem more likely (at 
higher temperature a new stem is much more likely to be 
removed) . 

The relationship between the probability distributions 
and the fixed-point attractor can be most clearly seen 
from the contour plot of log(/ bulk (i, j)) in Figure |5|b. It 
clearly shows the important roles played by / m i n and j 
in constraining the thickness of the next layer. When 
the previous layer is larger than /** the average stem 
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length must lie between Z m i n and j causing the crystal 
to become thinner. As j approaches l min the probabil- 
ity distributions become narrower and the difference in 
probability between the stem length being greater than 
j and smaller than j decreases until the fixed point I** 
is reached where the average stem length is the same as 



the thickness of the previous layer. This explains why I** 
is slightly larger than / m ; n . For j < I** the asymmetry 
of the probability distribution is reversed. These results 
confirm that the mechanism for thickness determination 
in the SG model isj/py similar to that which we found 
in the ULH modelJifM 




FIG. 5. (a) Probability distributions of the stem length in layer n given that the layer n+1 is of thickness j at T = 0.95T m . 
The labels give the value of j and the data points are for i = j. (b) shows a two-dimensional contour plot of log(/4(i, j)) at 
T — 0.95T" m . (The contours have a spacing of 0.5 and are measured with respect to the maximum in the probability. For 
clarity we only include the first 30 contours below the probability maximum; the probability does actually continue to fall 
rapidly in the blank bottom right corner.) The lefthand figures are for the outermost layer and the righthand figures for a 
layer in the bulk. The thicker line in the righthand panel of (b) displays the fixed-point attractor, Zbuik(j). 



We have also shown our results for f[(i,j) in Figure 
§. The main difference from f^ uik results from the fact 
that stems shorter than Z m j n can be present in the outer- 
most layer. The lack of this constraint causes the outer 
layer to be thinner than the bulk, and therefore the crys- 
tal profile to be rounded. The two modes of exponential 



decay at larger stem length that we noted for /£ ulk 
clearly also present in /{. 



are 



Another finding from our work on the ULH model 
was that the rate of growth actually slowed down as-tha 
thickness of a crystal converged to I** from aboveliSEij 
To determine whether this behaviour also holds for the 
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SG model we performed kinetic Monte Carlo simulations 
from initial crystals which had different thicknesses. As 
can be seen from Figure ^ the initial growth rate de- 
pends on the thickness of the initial crystal, but the final 
growth rate is independent of this initial thickness (on 
the right-hand side of the figure the lines are all parallel) 
because the thicknesses of all the crystals have now con- 
verged to I**. In agreement with our previous results the 
initial growth rate increases with the initial crystal thick- 
ness. Interestingly for a crystal with an initial thickness 
of 7 (less than Z m ; n ) the crystal initially melts causing 
the growth face to go backwards. Only once a fluctua- 
tion causes the nucleation of some thicker layers at the 
growth front can the crystal begin to grow again. The 
thinner the crystal the more difficult is this nucleation 
event. 



description to the picture presented above or whether the 
two pictures are incompatible in some respects. This re- 
examination is complicated by the fact that two slightly 
different arguments were presented. 

In the first the growth rate was decomposed into the 
termsEj 

G(i)=k+C 1 (i)-k-(l,i)Pi(l ) i), (19) 

which satisfy G = ^2 li G{i), The function G(i) has a 
maximum close to Zbuik (Figure ^) . It is expected that 

r bu ik = * G «/ G - (2°) 

i 

Sadler and Gilmer argue that the maximum in G(i) 
causes the crystal thickness to take the value Ztmik- 




FIG. 6. Growth from initial crystals of different thick- 
ness at T/Tm = 0.95. The initial position of the growth 
front is defined to be zero. The results were obtained from 
kinetic Monte Carlo simulations. The lines are labelled by 
the thickness of the initial crystal. 

Although we have only presented results for the two- 
dimensional version of the Sadler-Gilmer model in this 
section, we should note that in the three-dimensional ver- 
sion of the model we also found that the thickness of a 
crystal converged to the same value whatever the starting 
configuration of the crystal. 



B. Sadler and Gilmer's entropic barrier 

The description of the mechanism of thickness selec- 
tion that we give above is in terms of the fixed-point at- 
tractor, Zbuik(j)- It makes no mention of the growth rate 
(or a maximum in that quantity) in contrast to the LH 
surface nucleation theory and Sadlep,aiid Gilmer's own 
interpretation of the current model.tL2rt3 This of course 
raises the question, "Which viewpoint is correct?" In 
this section we examine Sadler and Gilmer's argument 
again to determine whether it provides a complementary 
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FIG. 7. Decomposition of the growth rate (solid line) 
into a barrier and a driving-force term using Equation (j2^) 
at T/T m = 0.95. 

Firstly, we should note that this explanation is very 
different from the maximum growth-rate criterion in the 
LH theory. In the LH theory the growth rate has a max- 
imum with respect to the thickness in a fictitious ensem- 
ble of crystals of different thickness, each of which grows 
with constant thickness. Rather, G(i) is a characteriza- 
tion of that one steady state where growth with constant 
average thickness occurs. 

Secondly as a characterization of the steady state 

Equation is correct. As G(i) is simply the rate of 
incorporation of stems of length i into thp-ytjjstal, at the 
steady state it should obey the equationiiirEil 



G(i) = GC hu]k (i). 



(21) 



Comparison of Figures |^b and ^ provides a simple visual 
confirmation of the validity of this equation. Equation 
( pp| ) simply follows from the above equation and Equa- 
tion ([l2] ) . The above equation also implies that explana- 
tions in terms of G(i) or C\, u ik(z) should be equivalent. 

Thirdly, although the steady state must be stable 
against fluctuations (the rate equations would not con- 
verge to it if it was not) G(i) does not in itself explain 
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how convergence to the steady state is achieved or how 
fluctuations to larger or smaller thickness decay away. 
G(i) only characterizes the steady state and not the sys- 
tem when it has been perturbed from the steady state. 
For example, the maximum in G(i) does not imply that 
crystals which are thicker than l^uik grow more slowly. 
Rather, as we showed earlier, these crystals initially grow 
more rapidly, but this growth causes the thickness to con- 
verge to I**. 

Sadler and Gilmer then continue by examining the 
form of G(i). Rearranging Equation ( |l9| ) leads toE^I 

This equation now has the form of a driving force term 
multiplied by a barrier term. The barrier term, C\{i), is 
the probability that the system is at the transition state 
and the driving-force term (fc + multiplied by the term 
in the brackets) is the rate at which this transition state 
is crossed. The two terms are shown in Figure ^ and it 
can be seen that their forms agree with the above inter- 
pretation. At small i, G(i) is small because of the small 
driving force and at large i the barrier term leads to the 
decrease in G(i). This leads to a maximum in the growth 
rate at an intermediate value of i when the barrier and 
driving force terms have intermediate values. 

It should be noted that the term we have labelled 
the driving force does not have exactly the behaviour 
one might [Expect. In contrast to Sadler and Gilmer's 
suggestion,Ej this term is not simply proportional to 
i — l-mm near to l m i n . For, as is evident from Equation 
(f2l|), G(i) can never become negative even at i < l m i n . 
This is because G(i) is not the rate of growth of a crystal 
of thickness i but the rate of incorporation of stems of 
length i when the system is at the steady state; there is 
always a small probability that a stem with length less 
than l m i n is incorporated into the crystal. 

The barrier term, C\, has a clear interpretation; it 
quantifies the effect of the rounding of the crystal growth 
front in inhibiting growth. The rate-determining step in 
the growth process is the fluctuation to a more unlikely 
squared-off configuration that has to occur before growth 
is likely to proceed. However, Sadler and Gilmer inter- 
pret C\ in terms of what they call an entropic barrier. 
Firstly, we note that a decay in C\(i) does not neces- 
sarily imply an increase in the free energy with i because 
C\ (i) is a steady-state probability distribution under con- 
ditions where the crystal is growing, whereas the Landau 
free energy for the outermost stem, Ai(i), is related to 
the equilibrium probability distribution through 

Ax(i) = A - fcriogC 1 oq (i), (23) 

where A is the Helmholtz free energy. The outermost 
layer may not have reached a state of thermodynamic 
equilibrium. 

Secondly, and more importantly, their explanation 
makes little mention of the limiting effect of the thickness 



of the previous layer. However, if we examine the contri- 
butions to Ci(i) from configurations where the thickness 
of layer 2 is j (Figure ||) it can be clearly be seen that 
these contributions decay sharply for i > j. This is be- 
cause it is unfavourable for a stem to extend beyond the 
previous layer. The combination of this decay in f{(i,j) 
(Figure |5|) with the rounding inherent in G^fJ) (Figure 
||b) causes C\(i) to decay as i increases. 
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FIG. 8. The solid line is Ci(i) and the dashed lines are 
the contributions to Ci from configurations where the thick- 
ness of layer 2 is j (Pi(i,j) — C2(j)f[(i,j)). Each dashed 
line is labelled by the value of j. T — 0.95T" m . . 

Sadler and Gilmer also considered a second approach 
to account for the thickness selection of lamellar crystals. 
In this approach a maximum in Gil), the growth rate of 
a crystal of thickness I, is sought. E£l However, by doing so 
this approach suffers the same weakness as the LH max- 
imum growth criterion. The implicit assumption that a 
crystal of any thickness continues to grow at that thick- 
ness is contradicted both by the behaviour .of .the current 
model (e.g. Figure ||d) and by experiment ESEll 

As with the previous approach the growth rate is di- 
vided into a barrier and a driving force term. The bar- 
rier term is C(Z), the probability that the system had a 
squared-off configuration (i.e. l n = I for all n). Sadler 
and Gilmer argue that C(l) is an exponentially decreas- 
ing function of I, because 'it seems likely the number of 
tapered configuijations will go up approximately expo- 
nentially with VIB 

Later, Sadler quantified this argument by calculating 
C(l) for the case that l n < Z ri+i :li3 

g(o-n 1+(t Vy < M > 

However, this result appears to be incorrect. We find: 

i-i 

C(l) = l[[l-(k-/k+y}. (25) 

3=1 
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The original and the present derivations are given in the 
Appendix £3 However, the two expressions do have a sim- 
ilar functional dependence on I and temperature, and so 
the error does not affect the substance of the argument. 
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FIG. 9. The probability of a squared-off configuration, 
C(l), at different temperatures as labelled. The diamonds 
mark the position of ibuik at that temperature. 

In Figure ||a we show the dependence of C(l) on I at a 
number of temperatures. It can be seen that C(l) does 
not decay exponentially, iuit after an initial decay reaches 
an asymptote at large Lt3 The origins of this behaviour 
in Equation (^5|) is easy to understand. For T < T m> 
k~ /k + < 1 and therefore 



lim 

l — >oo 



P(l) 



lim 

l—>oo 



1 



(k~/k J 



1. 



(26) 



The rate of decay to this asymptote depends on the value 
of k~ /k + . As the temperature is decreased, the mag- 
nitude of AF increases and therefore k~ /k + decreases 
(Equation [?]) and C(l) decays to its asymptote more 
rapidly. The curves in Figure y show just this behaviour. 

The physical significance of these results is clear. For 
thick crystals a squared-off configuration does not be- 
come any more unlikely as the thickness is increased. In 
particular, the free energy barrier to growth caused by 
rounding of the crystal profile does not grow exponen- 
tially with I as Sadler and Gilmer have argued. 

It should also be stressed that the rounding of the crys- 
tal growth front that is expressed by C(l) is not just a 
result of entropy. The number of configurations with an 
outermost stem of length i is independent of i. Round- 
ing of the crystal profile occurs because those fluctua- 
tions that lead to a stem near the edge of the crystal 
being thicker than the previous layer are energetically 
disfavoured compared to those which lead to a stem that 
is shorter than the previous layer. Again we see the cru- 
cial role played by the free energetic penalty for a stem 
extending beyond the previous layer. Sadler implicitly 
included this effect in the calculation of C(l) by impos- 
ing the condition l n < l n +\. Rounding of the crystal 



profile results from an interplay of energetic and entropic 
contributions to the free energy. 
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slope = -2(1 -T/T m ) 
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FIG. 10. (a) A x (i\j) for j = 10 and Ai(i). (b) Com- 
parison of the steady-state (dashed lines) and equilibrium 
(solid lines) values of the probabilities f[(i,j) and Ci(i). 
The equilibrium values were derived from Equations ( |29] ) 
and (|3C|). For f[{i,j) the lines are labelled by the value of j. 
T = 0.95T m 



C. The free energy of the outermost stem 

As the free energy function that is experienced by the 
outermost stem plays such an important role in both 
the LH and SG mechanisms of thickness determination 
here we calculate this quantity for the current model. 
This also allows us to examine the magnitude of non- 
equilibrium effects. 

Firstly, A\(i\j), the free energy of the outermost stem 
when the thickness of layer 2 is j, is simply 

A 1 {i\j)/e = 2iT/T m -2i+l i<j (27) 
A 1 (i\j)/e = 2iT/T m -i-j + l i>j. (28) 

For i < j the free energy decreases with increasing i be- 
cause the energy gained from crystallization is greater 
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than the loss in entropy. However, for i > j the new 
polymer units added to the crystal no longer gain any 
energy from interactions with the previous layer causing 
the free energy to increase. The example free energy pro- 
file, A\ (i\ 10), shown in Figure [To| a exhibits these features. 

From Ai(i\j) it is easy to calculate the value of f[{i,j) 
if the outermost stem is at equilibrium: 



exp(-Mi\3)/kT) 



(29) 



In Figure |10| b we show f[ eq {i,j) at a number of values 
of j. / 1 cq (i,j) is always peaked at j the position of the 
free energy minimum in Ai(i\j) and the asymmetry in the 
probability distribution (the stem length is more likely to 
be less than j than greater than j) reflects the asymmetry 
in the free energy profile. Comparison of the equilibrium 
and steady-state values of f{(i,j) shows that the two are 
in reasonable agreement for small j but the discrepancy 
increases with j. The steady-state values of f{(i,j) are 
always larger at small values of j because of the bias 
introduced by the successful initiation of a new stem be- 
fore equilibrium is reached for the outermost stem; this 
kinetic effect was noted in our earlier discussion of the 
form of f'{i,j). As j increases it becomes increasingly 
likely that a new stem is initiated before the current stem 
has reached a length j. 

We can also obtain C^ q (i) from fi eq (i,j)- 



(30) 



Using the steady-state values of C2 in the above equation 
we found that C^ q and C\ are virtually identical (Figure 
p^b) reflecting the good agreement between f[(i,j) and 
/1 q («, j) for small j. 

As C x q (i) is also related to A\(i), the free energy of a 
stem in the outermost layer of length i, through Equation 
@, it follows that 



Alii) 



kT log 

j 

/cTlogE 



C 2 (j)eM-Aim/kT) 
C2U) 

Ei< eM-Mm/kTY 



(31) 



where we have chosen the zero of Ai(i) so that j4i(1) = 
A x {l\j) = 2T/T m - 1. Figure [l(]a shows that A x {i) rises 
as i increases. There is a free energy barrier that increases 
with stem length. As with C\ this increase results from 
the combination of the rounding present in C2 and the 
free energy barrier to extension of a stem beyond the edge 
of the previous layer. 



are a number of differences in the behaviour of the mod- 
els. Firstly, the crystal profile in the SG model is always 
rounded whereas in the ULH model rounding does not 
occur at larger supercoolings. Secondly, /buik(j) in the 
SG model is a monotonically increasing function of tem- 
perature, whereas it increases at lower temperatures in 
the ULH model. This increase is because the free energy 
barrier for the formation of a fold (and hence the initia- 
tion of a new stem) becomes increasingly difficult to scale, 
so on average a stem continues to grow for longer. This 
effect is probably also partly responsible for the lack of 
rounding at larger supercoolings. One possible source of 
some of the differences between the models is the neglect 
of the energetic cost of forming a fold. In this section we 
consider the effect of non-zero e/ by presenting results 
for €f = e. 
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D. Effect of e 
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Although we found a similar mechanism ofpfciiickness 
selection in the SG model as in the ULH modcl,t§ll3 there 



FIG. 11. The dependence of the thickness of a layer in 
the bulk of the crystal on the thickness of the previous layer, 
j, and the temperature when e/ — e. In (a) the lines are 
labelled by T/T m and in (b) the lines are labelled by the 
value of j. 
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One effect of a non-zero ef is to increase l m [ n (Equa- 
tion |l3|) and so Tbuik is always larger than for ej = 
(Figure ||a). Furthermore, Tbuik is now a non-monotonic 
function of temperature; it gradually begins to increase 
below T = 0.8T m . We will explain the cause of this in- 
crease later in this section. The growth rate is slower 
than for e/ = 0, because the energetic cost of forming a 
fold makes initiation of a new stem more difficult. How- 
ever, the temperature dependence of G is fairly similar 
(Figure |b). 

The convergence of the thickness to its steady-state 
value is still described by a fixed-point attractor (Figure 
|IT|a) but the slope of the line is closer to 1 and so the 
convergence is less rapid. Furthermore, the temperature 
dependence of ^buik(j) is different (Figure |TT|b) . ?bulk(i) 
is now approximately constant at a value near to j over a 
wide range of temperature, whereas for ef = 0, ^buik(j) is 
approximately independent of?' at low temperature with 
a value close to i m i„ (Figure This difference is due 
to the greater difficulty of initiating a new stem when 
ef = e; the length of the stem is therefore much more 
likely to reach j. Similarly, comparison of Figure O with 
Figure || shows that the kinetic decay of /buik^'i)* 11 tne 
range l mnl < i < j is much reduced. As a result the 
steady-state and equilibrium probabilities for the outer 
layer are in closer agreement than for Cf = (Figure |l2|c) 
and the deviation of the rounding of the crystal profile 
from the theoretical prediction is reduced at low temper- 
ature (Figure ||c). 

The slopes of the free energy profile Ai(i\j) for i < j 
and i > j (Figure |io|a) have the same magnitude when 
T = 0.75T m . Therefore, below this temperature the 
probability distribution /( cq (i, j) is asymmetrical about j 
with i > j more likely. This asymmetry can also be seen 
for /b u jk(ijj) in Figure |l2|a. In the absence of the low i 
shoulder in /buik(*'i)> which is caused by the initiation 
of new stems before equilibrium in the outermost layer 
is reached, this asymmetry would lead to the divergence 
of the thickness. Instead, it just causes ^buik to rise mod- 
estly. By contrast, below T < 0.5T m there is no longer 
any free energy barrier to the extension of a stem beyond 
the previous layer and Zbuik increases rapidly. These ef- 
fects are not seen at e/ =0 because of the more rapid 
kinetic decay of f'(i,j) as i increases. 

Although using a non-zero ef makes the behaviour of 
the SG and ULH models in some ways more similar, dif- 
ferences remain. For example, the crystal profile is still 
always rounded (Figure ||c). The differences probably 
stem from the differences in the number of growth sites 
in the two models. In the ULH model the crystalline 
polymer is explicitly represented and there is only two 
sites in any layer at which growth can occur — at either 
end of the crystalline portion of the polymer. In contrast 
growth can occur at any stem position at the growth 
front in the SG model. Consequently, in the ULH model 
most of the stems in the outermost layer must be longer 
than Z m i n if the new crystalline layer is to be stable, thus 
making the rounding of the crystal profile less than for 



the SG model. 
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FIG. 12. Probability distributions of the stem length in a 
layer in the bulk of the crystal given that the previous layer 
is of thickness j at (a) T = 0.6T m , and (b) T = 0.95T m . In 
(a) and (b) the labels give the value of j and the data points 
are for i = j. (c) Comparison of the steady-state (dashed 
lines) and equilibrium (solid lines) values of the probabilities 
Ci(i) and f! The equilibrium values were derived from 

Equations ( p9[ ) and (|30[). For f'(i,j) the lines are labelled 
by the value of j. T — 0.95T m . All results are for tf = e. 
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We can also reach this conclusion if we note that 
in some ways the ULH model is similar to the two- 
dimensional SG model but where the slice represents 
a polymer growing along the surface rather than a cut 
through the crystal perpendicular to the growth face.E£l 
For this variant of the SG model, Figure |5ja would rep- 
resent the dependence of the average stem length in the 
outermost layer on the distance from the growing end of 
the crystalline portion of the polymer. Only near to this 
growing end could the stem length be less than Z m ; n . 



IV. CONCLUSION 

In this paper we have reexamined the mechanism of 
crystal thickness selection in the Sadler-Gilmer (SG) 
model of polymer crystallization. We find that, as-.with 
the ULH model that we investigated previously,li§li3 a 
fixed-point attractor describes the dynamical conver- 
gence of the thickness to a value just larger than l m i n 
as the crystal grows. This convergence arises from the 
combined effect of two constraints on the length of stems 
in a layer: it is unfavourable for a stem to be shorter 
than Z m j n and for a stem to overhang the edge of the pre- 
vious layer. There is also a kinetic factor which affects 
the stem length: whenever a stem is longer than / m ; n the 
growth of the stem can be terminated by the successful 
initiation of a new stem. This factor plays an important 
role in constraining the stem length at larger supercool- 
ings when the free energy of crystallization is larger and 
the barrier to extension of a stem beyond the previous 
layer is therefore smaller. 

This explanation of thickness selection differs from that 
given by Sadler and Gilmer. They realized that in or- 
der to explain the temperature dependence of the growth 
rate in their model they required a 'barrier' term that irt. 
creased with the crystal thickness at the steady state.El 
They correctly identified that the rounding of the crystal 
profile has just such an inhibiting effect on the growth 
rate. However, they then also tried to explain the thick- 
ness selection in terms of the growth rate in a manner 
somewhat similaip-to the maximum growth rate approach 
in the LH thcory£3 This approach has a number of prob- 
lems. Firstly, G(i) (Equation |is| ) merely characterizes 
the steady state. A maximum in that quantity does not 
explain how or why convergence to that steady state is 
achieved. G(i) is just the rate of incorporation of stems 
of length i into the crystal at the steady state. Equa- 
tion ( |2l| ) trivially shows that it must have a maximum 
at the most common stem length in the crystal. By con- 
trast, Sadler's alternative maximum growth argument in 
terms of C{1) has potentially more explanatory power, 
but it falls into the LH trap of incorrectly assuming that 
a crystal of any thickness can continue to grow at that 
thickness. 

Secondly, it is somewhat misleading to label the effect 
of the rounding of the crystal profile as due to an en- 



tropic barrier because it is due to both the energetic and 
cntropic contributions to the free energy of the possible 
configurations. In particular, it ignores the important 
role of the free energy barrier for the extension of a stem 
beyond the edge of the previous layer. 

We therefore argue that our interpretation in terms 
of a fixed-point attractor provides a more coherent de- 
scription of the thickness selection mechanism in the 
SG model. The added advantage of this pictures is its 
generality — we have shown it holds for both the SG and 
ULH models. Therefore, the mechanism is not reliant 
upon the specific features of these models. 

There is also experimental evidence that a fixed-point 
attractor describes the mechanism of thickness selection 
in lamellar polymer crystals. The steps on the lamel- 
lae that result, torn a change in temperature during 
crystallizatior£j'c3 show how the thickness converges to 
the new /**. The characterization of the profiles of these 
steps by atomic-force microscopy could enable the fixed- 
point attractors that underlie this behaviour to be ex- 
perimentally determined. As well as testing our theo- 
retical predictions, such results would provide important 
data that could help provide a firmer understanding of 
the microscopic processes involved in polymer crystalliza- 
tion. For example, although the mechanism of thickness 
selection in both the ULH model and the SG model is 
described by a fixed-point attractor, the form of the tem- 
perature dependence of ?bulk(i) has significant variations 
depending on the model used and the parameters of the 
model (compare for example, Figures 
paper and Figures 6a and 13a in Ref. 



and [Tl] b in this 



ACKNOWLEDGMENTS 

The work of the FOM Institute is part of the research 
programme of 'Stichting Fundamenteel Onderzoek der 
Materie' (FOM) and is supported by NWO ('Nederlandse 
Organisatie voor Wetenschappelijk Onderzoek'). JPKD 
acknowledges the financial support provided by the Com- 
putational Materials Science program of the NWO and by 
Emmanuel College, Cambridge. We would like to thank 
James Poison and Wim de Jeu for a critical reading of 
the manuscript. 



APPENDIX A: CALCULATING C(L) 

In Ref. |lj Sadler presented an expression for the prob- 
ability that a crystal of thickness I has a squared-off pro- 
file (i.e. l n = I for all n). To do so he assumed that 
only tapered configurations were possible, i.e. l n < l n +i, 
and that the system is at equilibrium. The calculation 
requires that the ratio of the Boltzmann weight of the 
squared-off configuration to the partition function, Z, be 
calculated. If we use the free energy of the squared-off 
configuration for our zero, C(l) = 1/Z. 
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FIG. 13. Schematic diagram of the crystal profile to show 
(a) Sadler's approach and (b) our approach for calculating 
the probability of a squared-off configuration. The shaded 
area is proportional area to the free energetic cost of a con- 
figuration. 

Sadler argues that C(l) can be calculated using a 
random-walk analogy as in Figure [l3|a. Each possible 
configuration is represented by by a walk of horizontal 
and vertical displacements from i — 1 to i = I. To achieve 
a squared-off configuration I — 1 successive vertical dis- 
placements from i = l have to be made. As a horizontal 
move would involve a loss in free energy of —(I — i)AF, 
where AF is the free energy of crystallization per poly- 
mer unit, the ratio of the probabilities of horizontal and 
vertical moves is 



Ph 

Pv 



exp((l-i)AF/kT), 



It simply follows from ph + p v = 1 that 



1 



Pv 



1 + exp((Z - i)AF/kT) ' 



(Al) 



(A2) 



Therefore, the probability of I — 1 successive vertical steps 
is 



i-i 



1 



cxp((7 - i)AF/kT) 



(A3) 



Substituting the exponential by a ratio of rate constants 
(Equations [7j and |9|) and using j = I — i gives 



l-i 



c « = n r 



3=1 



(k-/k 



(A4) 



However, the method of calculation is flawed. A proper 
calculation of the partition function requires that all pos- 
sible configurations be considered. We can characterize a 



general configuration by the set of integers {j±, . . . , ji-i} 
(Figure |l3|b). The free energy of a configuration with re- 
spect to the squared-off configuration can be calculated 
by summing the shaded rectangles in Figure 13 d; this 
gives — ^2\=iji^ ~ i)AF. The partition function can 
then be written as 



'2-1 



Z 



exp[J2ji(l-i)&F/kT 

jl, ■■■,3l-l=0 



(A5) 



As an exponential of a sum is a product of exponentials 

2—1 oo 

Z = l[Y, c MMl-i)&F/kT). (A6) 

i=l j«=0 

The sum is a geometric series so 



«-n r 

i=l 



1 



exp((l -i)AF/kT)' 



(A7) 



Using C(l) — 1/Z and again substituting for the expo- 
nent and for i gives 



2-1 



c(o = n [i (k-/k + y] 

i=i 



(A8) 



In this approach we can also calculate a number of 
other quantities. For example, C\{i) — Z±(i)/Z, where 
Z±(i), the contribution to the partition function from 
configurations with an outer stem of length i, is given 
by 



2-; 



Zl ( l ) = (k-/k+y-*i[ T 



3 = 1 



{k-/k+y 



(A9) 



for i < I - l and by Z x (l - 1) = (1 - k~/k + ) 1 and 
Z\iX) = 1. l\ can then be calculated using Equation 
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